Black Holes and Wormholes in spinor polariton condensates 
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We propose a new system for the study of event horizons and black holes - a Bose-Einstein 
condensate of exciton-polaritons. Hawking radiation is observed in numerical experiment. A closed 
horizon is obtained in 2D. We simulate inter- Universe and intra- Universe wormholes capitalizing 
on both the spinor nature of polariton condensates and the spin dependence of polariton-polariton 
interactions. 
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It happens quite often in Physics that striking simi- 
larities are found between systems, which from the first 
glance have absolutely nothing in common. Sometimes, 
such similarities can be exploited to perform laboratory 
studies on accessible objects similar to inaccessible ones. 
One of the most recent examples is that of Klein tunnel- 
ing in graphene[l|: the quasi-particles in a solid-state ob- 
ject obey the same mathematical equations as very high- 
energy relativist ic particles. Astrophysics allows even less 
laboratory studies than high-energy physics: scientists 
are restrained to the objects in the Universe proposed by 
the Nature, and these are studied from very far. Thus, 
having a desktop version of a black hole (BH) would be 
even more useful. Once again, the Physics of the small 
comes to the aid of the Physics of the large. The anal- 
ogy between the equations describing the excitations of a 
Bose-Einstein condensate (BEC - a peculiar state of mat- 
ter exhibiting quantum properties at macroscopic scales) 
and the metrics of the curved space-time has been no- 
ticed about a decade ago[2|. Since then, the scientists 
have managed to experimentally observe the event hori- 
zons in atomic BECs[3|. 

However, such atomic condensates are still a bit far 
from being a convenient laboratory tool, because they 
require ultra-low temperatures for their formation, and 
the measurements of the distributions inside the conden- 
sates are relatively complex to carry out. Here the solid- 
state physics comes into play with exciton-polaritons Q: 
hybrid quasi-particles that share the properties of both 
light and matter, formed in microcavities in the strong 
coupling regime. On one hand, their excitonic part allows 
them to interact both with themselves and the surround- 
ing phonon bath, providing efficient relaxation and for- 
mation of a BEC [5|- On the other hand, their photonic 
part attributes them a very small effective mass which 
allows BEC even at room temperature The finite 
(and short) lifetime of polaritons turns into an impor- 
tant advantage, simplifying all measurements, because 
the decay of the condensate means the emission of pho- 
tons from the cavity, and the distribution of emitted pho- 
tons gives direct information on the polaritons distribu- 
tion function, on their dispersion and spatial evolution, 



that is, the whole wavefunction ip (r, t) up to a constant 
global phase (using interferometry ) . Finally, the spin 
structure of polaritons (two allowed spin projections 
on the growth axis) implies possible new effects due to 
the vectorial nature of the condensates. These partic- 
ularities favor polaritons with respect to other systems 
proposed for the simulation of BHs, including the optical 
ones based on metamaterials [8l4l0|. 

In this paper, we propose polariton condensates as a 
convenient tool for the study of sonic BHs, event hori- 
zons and Hawking radiation. We begin with a simpler 
case of a ID scalar condensate and then consider a 2D 
system with a closed horizon. Finally, capitalizing on 
spin, we demonstrate different configurations containing 
wormholes, including one allowing the signal to be trans- 
mitted faster than the speed of sound. 

Event horizon in ID flow It has been understood quite 
a long time ago in hydrodynamics, that an event horizon 
can appear if the flow speed increases and becomes larger 
than the speed of sound [ll|. Indeed, the excitations in 
the flowing medium propagate with the speed of sound 
(in the linear approximation, that is, long-wavelength 
limit), and therefore in the laboratory frame they are 
unable to go against the flow if its speed is too high. A 
great research effort in this domain has recently culmi- 
nated with the observation of stimulated Hawking emis- 
sion in water [l2| . 

The weak excitations in the BECs (called Bogolons) 
can also exhibit linear dispersion in the long- wavelength 
limit, with the speed of sound given by c s — ^an/m, 
where a is the interaction constant (in the Born approx- 
imation), n is the density and m is the mass of bosons 
composing the condensate. The finite lifetime of polari- 
tons provides a natural way to vary the speed of sound. 
A propagating condensate can be injected locally by res- 
onant or non-resonant pumping [13[ . Its density is bound 
to decrease with the distance from the pumping spot, and 
the speed of sound will therefore decrease as well, whereas 
the propagation speed will remain constant or may even 
increase, if the condensate is accelerated by some poten- 
tial ramp (for example, by its own self- interactions). At 
some point the two speeds become equal, defining the 
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position of the event horizon of the BH. 

To describe accurately the polar iton condensate, we 
use a set of standard equations for the photonic 
tjjph(x, t) and excitonic ij) ex (x, y, t) mean fields coupled 
via the light matter interaction Qr = 15 meV, account- 
ing for the two allowed spin projections a = ±1: 
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Here m p h = 3.6 x 10 _5 rao, m ex = OAtuq and mo are 
the cavity photon, the quantum well exciton and the free 
electron masses respectively. The lifetime of the parti- 
cles are r p h = 25 ps and r ex = 300 ps. U(x) is the po- 
tential (e.g. disorder). H e ff(x,y) accounts for the pos- 
sible effective magnetic field coupling the two spin com- 
ponents. The interaction constants a\ = 0.02 meV//mi 2 
and \ct2 1 << \ol\ | [H| correspond to parallel and anti- 
parallel spins. Finally, P CJ (x 1 y 1 t) is a photonic pump 
tuned in quasi-resonance with the lower polar iton branch. 

First, we consider a ID case neglecting the spin de- 
gree of freedom. The polariton flow is resonantly and lo- 
cally injected by a pumping laser located close to x = 
on the figure QJa). The medium is assumed to exhibit 
small structural disorder mainly due to the etching of 
the ID wire cavities. It is modeled by a random series 
of delta-peaks separated by 1 fim on average. All re- 
sults shown in FigJT]are averaged over 100 disorder real- 
izations. The FigQJa) shows the polariton density n(x), 
the visible decay is mainly due to the finite lifetime of po- 
laritons. The speed of the sound c(x) decreases together 
with the density. On contrary, the speed of the flow v(x) 
is increasing because of the self-interactions within the 
condensate. In the left part, the flow is subsonic and 
cannot be scattered by the disorder (superfluid). In the 
right part, the flow is supersonic. The exponential decay 
is induced not only by the life time but also because of 
the Anderson like-localization in the disorder. The two 
regions are separated by a horizon at v = c. In the su- 
personic region no excitation can propagate towards the 
horizon. The generation of Hawking emission on the hori- 
zon is demonstrated in the FigOJb). Indeed as recently 
proposed in Refs.[l5|, [Hj, emission of Hawking phonons 
means correlated density perturbations propagating on 
both sides of the horizon. Hawking emission can there- 
fore be detected using the density-density correlation ma- 
trix (x,x f ) = (n (x) n (x')) / {n (x)) (n(x')). This ma- 
trix, averaged on disorder, is shown at t = 200 ps in 
FigOJb). Indeed, as expected, the "Hawking tongues", 
indicating positive correlations, are extending from the 
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FIG. 1. (Color online) a) Propagation of a polariton conden- 
sate in ID: density (black solid line), flow speed (red dashed 
line), sound speed (green dash-dotted line). The horizon is in- 
dicated by a blue vertical line, b) Density-density correlation 
matrix g( 2 \x,x') at t — 200 ps. Red lines are a guide for the 
eyes, indicating the Hawking tongues (positive correlations) 
extending from the main diagonal. 



horizon position, marked by the red dashed lines. We 
underline that we do not need to introduce quantum fluc- 
tuations to seed Hawking emission thanks to the presence 
of the disorder potential and to the finite lifetime which 
broaden the states in momentum and frequency. 

2D black holes The finite lifetime of polaritons allows 
to organize persistent flows, as shown in ref. [l3|, Hl|, [22[ 
for the cases of quasi-resonant or non-resonant pumping. 
This particular property makes the formation of a closed 
event horizon in 2D possible. This is much more com- 
plicated with atomic condensates, where only ID config- 
urations have been considered With polaritons, one 
needs to pump a region around a large-scale defect in 
the microcavity mirrors possessing a lower quality fac- 
tor. Polariton-polariton interactions will create a per- 
sistent flow converging into the defect region, where the 
density is always lower due to the shorter lifetime. 

Figure [2] shows the results of a realistic 2D simulation 
with pulsed spatially homogeneous pumping. The pho- 
ton density \^ p h\ 2 at the time t = 8 ps is plotted as a 
function of coordinates (panel a) . The defect region with 
a shorter lifetime is located at the origin ((#, y) = (0, 0)), 
while the disorder is neglected. The density inside the 
defect region decreases faster than outside, and the re- 
pulsive interactions make polaritons propagate towards 
the center of the figure. The event horizon at that time 
is marked with a dashed line. The Hawking radiation is 
seeded by the non-equilibrium spatial distribution. It can 
be observed on the panel b) as density waves propagating 
inwards inside and outwards outside the horizon. 

Hawking radiation in BECs has been largely studied 
in previous^ works (see e.g. Q for low- wavelength limit 
and M ' 



24j for the general case) and is not the main 
object of the present paper, which is why we restrain 
ourselves to the numerical demonstration of this effect, 
even though the possibility of experimental study of the 
Hawking radiation was one of the main reasons for the 
beginning of the activity on acoustic horizons [11 1. 
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FIG. 2. (Color online) a) 2D BH around a defect (dark region) 
in a polariton condensate. Dashed line shows the event hori- 
zon, and red arrows point the directions of the propagation of 
the Hawking radiation on both sides of the horizon, b) Wa- 
terfall density plot at different times. Propagating Hawking 
phonons are marked by the red dotted lines. 



Wormholes If a single scalar condensate is a model of 
a Universe which might contain some BHs, it is natural 
to map two spinor components to 2 different Universes. 
They can be completely decoupled from each other if 
there are no interactions between the particles of different 
spins. Adding a magnetic field can provide a coupling 
between these two Universes, making possible simulation 
of wormholes 



In astrophysics, inter-Universe wormholes are the pairs 
of singularities located in different Universes and con- 
nected together. Using such wormhole, one can pass 
from one Universe to the other. A more interesting situa- 
tion, when both holes connected together are in the same 
Universe, is called an intra-Universe wormhole. Such 
wormhole can connect two distant regions of space with a 
tunnel much shorter than the distance between the two, 
which might allow faster-than-light travel [l8| . 

We start with the simpler case of an inter-Universe 
wormhole. We consider a ID quantum wire, as in [l3| . 
The idea is to first create a closed BH bordered by two 
event horizons and to connect the latter with a white hole 
in the other spin component using a local effective mag- 
netic field H e ff, induced by the energy splitting exist- 
ing between linearly polarized (TE or TM) eigen modes. 
H e f f be controlled by varying the width of the wire [l9[ , 
or by applying an electric field [20| . 

In general relativity, a key concept is the propagation 
of signals, whose speed can never exceed that of light 
in the vacuum. The propagation of phonon wavepack- 
ets across the event horizons in BECs has already been 
studied, for example in refs. TBI, [26|. In our model sys- 
tem, the "signals" will be grey solitons [25| created in a 



single spin component [27|, [28| (half-solitons) of the po- 



lariton condensate (e.g. <r + , representing our Universe). 
Such solitons have a lot of properties similar to those 
of relativistic particles [25[ except that their mass 
is negative (because they are actually holes rather than 
particles): m so i = mo/ yl — v 2 / 'c 2 ; their size is given 
by I sol = £/\A — v 2 /c 2 wnere v is their velocity, c s is 
the speed of sound and £ is the healing length of the 



condensate. For attractive interactions, the mass can 
of course be positive. We create pairs of half-solitons 
thanks to a short time-dependent pulsed potential act- 
ing on a single spin component U(x,t) = L r oexp(— (x — 
xq) 2 /w 2 ) exp(— (t — to) 2 /r 2 ). As the speed of a soliton is 
related to its depth by v = c s n(0)/n where n(0) is the 
density at its center and n in its surrounding [25], we need 
a weak enough potential U(x,t) to excite shallow solitons 
that will be able to travel at speeds close to c s in both 
spin components (marked as ci,C2 below). 

Figure [3] shows the scheme of a numerical experiment 
with a single wormhole (panel (a)) and the results of 
the simulations (panels b,c). A ID polariton wire is 
cw-pumped by two spatially separated, quasi-resonant, 
<j + polarized lasers allowing the formation of a steady 
state flow and specific density distribution sketched on 
the panel (a). Polaritons flow away from the pumping 
regions and a closed BH is formed at the crossing of the 
flows. An effective magnetic field converting <j + to <j~ is 
present in the BH region. The a~ density therefore shows 
a maximum expelling excitations outside from the central 
region which corresponds to the formation of a white hole 
in the a~ universe. After the steady state is obtained, a 
pulsed potential is applied in the <j + -component (panel 
b) at t = 5 ps at the left of the left horizon, creating 
a pair of propagating half-solitons. One of these prop- 
agates freely to the left, whereas the other enters the 
BH and remains partly guided inside. The effective mag- 
netic field converts a part of this soliton into the other 
spin component (a - , panel c). The soliton is then able 
to cross the horizon of the white hole, propagating away 
together with the flowing condensate. One can also see 
that short- wavelength perturbations are still able to cross 
the horizon of the BH in any direction. The holes pos- 
sess an internal structure, which we do not discuss here 
in details. The speed of the flow is zero at the center 
and the BH is in fact composed by two narrower BHs 
surrounding a subsonic region. 

In the last part, we propose a scheme for an intra- 
Universe wormhole allowing for the transfer of a half- 
soliton with an apparent velocity faster than the speed 
of sound of its original universe. Such intra-Universe 
wormhole is based on two inter-Universe wormholes sim- 
ilar to the the ones previously described, but connecting 
the Universes in opposite directions. The scheme (fig 
4a) shows the proposed density profiles with the two 
wormholes. Dashed line indicates the propagation of a 
half-soliton. A <j + half-soliton is generated at the left. It 
enters in the <j + BH where it is converted in a a~ soliton, 
which is ejected by the white hole part of the wormhole. 
Then, it travels in the a~ component between the two 
wormholes with a velocity close to C2- It then reaches the 
second wormhole, which is a BH in the cr~ component, 
where it is captured and converted to a <j + soliton ejected 
from the white hole. The average velocity of this soliton 
is close to C2 which can be larger than ci, the speed of 
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FIG. 3. (Color online) a) Scheme of a wormhole between the 
cr + (red) and cr + (blue) Universes. Arrows show the directions 
of the flow in the two components. Dashed lines mark the 
event horizons in both components, b) Results of numerical 
simulation: \ip(x,t)\ 2 for both spin components with black 
dashed lines showing the boundaries of the black/white holes 
and white arrows indicating the propagation of signals. 



sound in the <r + universe. The results of corresponding 
simulations are presented in panels (b) and (c). 

The pumping is cw quasi-resonant with inhomogeneous 
elliptical polarization, providing the density profiles close 
to fig. 4(a). A pair of half-solitons is created at t = 5 ps 
in the <j + -component at x = (panel b). The "reference" 
half-soliton propagates to the left with the speed limited 
by ci and arrives to the edge at around t = 30 ps. The 
half-soliton falling inside the BH converts into the o~~ 
component, gets out of the white hole in a~ (panel c) 
and propagates with a higher speed, limited by c 2 > c\. 
This half-soliton arrives to the second wormhole and is 
converted back into the <j + -component, appearing there 
at t = 25 ps (marked by the green circle). The two 
events are marked with black horizontal lines, and the 
time difference between them is At ~ 5 ps. 

Conclusions Spinor polariton condensates, being rel- 
atively easy to produce and manipulate, can be used 
for the simulation of astrophysical objects, such as black 
holes and wormholes. A 2D black hole with a closed event 
horizon can be simulated. Effective magnetic fields, well 
known as the cause of non-trivial spin dynamics of polari- 
tons, can be used to organize the coupling between the 
black holes and white holes in the two spin components. 
A system of two wormholes in opposite directions allows 
one to organize "faster-than-sound" signal propagation. 
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showing density profiles; b) cr + density; c)a~ density. Dashed 
horizontal lines mark the arrival of the two signals (green 
dashed circles), with At the time difference between them. 
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